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It is argued that the pairing symmetry realized in a ferro- 
magnetic superconductor UGe2 must be a non-unitary triplet 
pairing. This particular state is free from the Pauli limitation 
and can survive under a huge internal molecular filed. To 
check our identification we examine its basic properties and 
several experiments are proposed. In particular, the external 
field is used to raise Tc by controlling the internal spontaneous 
dipole field. 

PACS numbers: 74.20.-z,75.50.Cc,74.70.Tx 

Ferromagnetism (FM) and superconductivity (SC) 
are thought to be basically mutually repulsive. Since 
Ginzburg [0 points out a possibility of its coexistence un- 
der the condition that the magnetization is less than the 
thermodynamic critical field, there are many experimen- 
tal investigations performed, starting by Matthias et al 
1^ who consider impurity ferromagnetism in a supercon- 
ductor. Although the coexistence between antiferromag- 
netism and superconductivity are rather easy to realize 
and actually observed in several compounds ^ because 
the antiferromagnetic moments spatially averaged over 
the SC coherent length vanish, the ferromagnetic case is 
difficult to realize. Rare exceptions of this case are rare- 
earth ternary compounds HoMogSg and ErRh4B4 where 
in a narrow temperature region just below the Curie tem- 
perature TpM the coexistence of FM and SC is attained 
Q . When the rare-earth 4f moments completely align at 
lower T, SC is wiped out by a strong internal field. So 
far there is no known SC which can fully sustain such a 
large molecular field. 

As for the theoretical developments, Anderson and 
Suhl 1^ discuss the cryptoferromagnetism as a possible 
coexistence phase where the FM is modified to a long- 
period modulated spin structure. This modified RKKY 
interaction is mediated by superconducting pairs. Simi- 
larly to this idea, Blount and Varma 1^ also propose a 
magnetic spiral phase coexisting with SC to gain electro- 
magnetic dipole interaction between localized moments. 

The recent discovered material UGe2 seems to be 
difficult to understand in terms of these theories so far 
proposed and seems to require a novel concept to inter- 
pret the data: Because previous theories assume that two 
groups of electrons are distinguishable and clearly sepa- 
rated spatially. FM comes from the well localized 4f elec- 
trons while SC pairs are formed by conduction electrons. 



Here the situation is much more intricate; we cannot sep- 
arate two groups in a well defined way from the outset, 
where 5f electrons from U-atoms play double roles both 
for FM and SC. It is a quite interesting problem how to 
describe its double roles played by 5f electrons microscop- 
ically, here we confine our discussions to a phenomeno- 
logical level without going into pairing mechanism, which 
is closely related to it. 

At ambient pressure UGe2 is an itinerant metallic fer- 
romagnet whose Curie temperature rFM=52K and the 
spontaneous moment ^ lAfiB /U-atom. The easy axis 
is the a axis in orthorhombic crystal |^. Upon increas- 
ing pressure P, both TpM and the spontaneous moment 
moderately decrease gradually and at P ~lGPa the SC 
starts to appear and the SC transition temperature in- 
creases, taking a maximum Tc ~ 0.8K above which FM 
already sets in at a high T of TpM ^ 30K whose mo- 
ment ^ l/^s/U-atom. Further increasing P, both TpM 
and Tc drop down to disappear almost simultaneously 
around P ~1.7GPa. Thus the SC region is entirely cov- 
ered by FM. The SC is confirmed to coexist with FM 
by neutron experiment This phenomenon is diffi- 
cult to understand in terms of the previous framework 
based in a conventional singlet pairing state. The in- 
ternal ferromagnetic molecular field coming from the or- 
dered magnetic system through the exchange interaction 
amounts to an order of few hundred tesla (T) in view of 
the exchange splitting ~70meV for the up spin and down 
spin Fermi surfaces . According to Tsutsui et al pd| ] 
the hyperfine field probed by U Mossbauer specroscopy 
is '^240T. This huge exchange field apparently excludes 
not only any singlet pairing category, but also certain 
forms of triplet pairing category, namely, unitary triplet 
states. These are all limited by the Pauli paramagnetic 
field Hp -1.3rc(T). Anderson and Suhl §, and Blount 
and Varma ^ theories are not applicable here, and Fulde 
and Ferrell (FFLO) state only slightly enhances Hc2 by 
at most 10% imi, but definitely not possible for >100T. 

In fact, only possible pairing symmetry under such a 
strong internal field is non-unitary triplet state, which is 
free from the Pauli limit |T^. Thus, we investigate this 
possibility in light of the present material UGe2 , examine 
the available data and predict some of the interesting 
phenomena associated with this non-unitarity |Q . 

A non-unitary triplet state is described by the order 
parameter A(k) = «(d(k) • (j)ay in a 2x2 matrix form. 
The three dimensional complex vector d(k) fully char- 
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acterizes the triplet pairing state. If d(k) is a complex 
number, the product A(k)A(k)'l' = |d(k)p(To + i(d(k) x 
d*(k)) • a is not a multiple of unit matrix and A(k) be- 
comes non-unitary. Thus in a non-unitary state time re- 
versal symmetry necessarily is broken spontaneously and 
spontaneous moment m(k) oc id(k) x d*(k) appears at 
each k point, yielding the macroscopic averaged moment 
< m(k) > provided that the its Fermi surface average 
non-vanishes. There are two kinds of non-unitary state 
among several possible forms, depending on its orbital 
structure; one is bipolar type and non-bipolar one. In 
the former (latter) the real and imaginary parts of the 
complex d(k) vector are ascribed by different (same) or- 
bital function(s). The average spontaneous moment in 
the latter (former) is nonvanishing (vanishing) in general 
by symmetry. Since the additional Zeeman magnetic en- 
ergy is gained in the non-bipolar, in the following we only 
consider the non-bipolar state. 

We derive an appropriate Ginzburg-Landau (GL) free 
energy density functional allowed in the presence of the 
FM described in terms of the complex order parame- 
ter 7/ of the coefficients of the d vector. The macro- 
scopic ferromagnetic order parameter M is proportional 
to the magnetization which is estimated as ~0.2T in the 
ambient pressure The complex order parameter 77 
is coupled linearly and quadratically with M, namely, 
Z7M ■ ff X ff* and Al'^ff ■ ff* respectively. Thus it is con- 
venient to write fj = {rj^, r/y, 0) or r]± = -^{rjx ± i^y) by 
choosing M = (0,0, M), implying that the Cooper pair 
spin orientation points to the M direction. We choose a 
coordinate system: x \\ b,y \\ c,z \\ a where the magnetic 
easy axis is the a axis. In view of the lattice constant a 
and c are almost same, only 2% difference, we regard the 
crystal structure as tetragonal where the (a, c) forms the 
basal plane. 

So far we have discussed the possible coupling terms 
with FM allowed by symmetry in order to explicitly see 
the effects of the FM formation. If we further consider 
the microscopic origin of SC, which appears under the 
Zeeman split Fermi surfaces (FS), leading to the distinct 
transition temperatures and the different fourth order 
terms for 'q± from the outset. Thus we finally arrive 
at the following generic GL form appropriate to UGe2 : 

hulk = + Oi-\V-? 

+ i/?+h+|4 + /3|r;+nr;_|2 + i/3_h_|4 (1) 

where we introduce a± = o^iT — T^) and /?± > for 
the stability condition of the system. In the following 
is identified to the observed Tc below which 77+ becomes 
non- vanishing. 

It is easy to derive the condition for the second SC 
transition to occur: a^T" / a'^T^ > f3/(3+ at Tc2 = 
{P+a-^T- - (3a^T+)/{(}+aQ - I3a^) at which the re- 
maining order parameter rj- starts to appear. Thus there 



is a chance to find the double transition near the critical 
pressure 1.7GPa where FM is suppressed and T~ will ap- 
proach r+. We note that this transition is analogous to 
the Ai-A phase transition in superfluid '^He described by 
Ambegaokar and Mermin . Near the critical pressure 
where M is small, we can estimate relative changes of the 
two transitions by expanding the density of states N^ep) 
at eF, giving rise to Tei,2 = r,o(l±M^^ In ^). The 
derivative of the densities of states N'^ep) and the energy 
cutoff UJ. 

The quasi-particle excitation spectrum for a non- 
bipolar type d(k) — f](f)(]i), where the spin ff and 
the orbital 0(k) parts are separable, has two branches: 
E^(k) = v/e2(k) -f A2(k) (cr = ±). The gap functions 
are given by A±(k) = |77j;±?7j^|(/)(k) = |77±|0(k). Thus one 
of the branches, say, A_(k) vanishes identically. That 
means that on the spin down Fermi surface there is no 
superconducting gap formed, leaving it normal. 

The nodal structure associated with non-unitarity 
leads to several observable predictions without specify- 
ing particular orbital function 0(k). The total density 
of states is given by N{E) = N+{E) + N^{E) where 
N„{E) = T,i^S{E- E„{k)). At the SC transition temper- 
ature, only the spin up FS opens the gap, thus the jump 
of the specific heat is substantially reduced from the BCS 
value (1.43). In the lowest T the specific heat becomes 
C{T) = 7_T where 7_ is the densities of states at the 
FS for the spin down band. According to an estimate by 
the spin polarized band structure calculation 7_ is 
substantial compared with the total 7 (=7+-|-7_). Thus, 
it is quite observable. The existence of the 7_ term gov- 
erns the lowest T thermodynamics, such as a T-linear 
term in thermalconductivity and the I — V characteris- 
tics of quasi-particle tunneling between a normal metal 
and UGe2, and etc. 

We notice that if the second transition Tc2 really takes 
place at lower T, there is no residual 7_ contribution to 
the thermodynamic quantities mentioned above since the 
remaining spin down FS is also gapped by 77- . 

As for the spin susceptibility Xi {i — a,b and c) probed 
by the Knight shift experiment, we naively expect that Xb 
and Xc (Xa) must decrease (unchanges) below Tc because 
the d(k) vector lies in the (6, c) plane, but remain a finite 
value, corresponding to 7_. Note that Xb and Xc may not 
decrease appreciably below because these transverse 
susceptibilities in FM are known to be determined by the 
whole bands and insensitive to the gap formation near 
FS. The direction of the d(k) vector is strongly locked 
to M//a where UGe2 is known to have strong easy axis 
type of the magnetic anisotropy and the magnetization 
is Ising-like Q]. 

Having analyzed the basic thermodynamic properties 
of non- unitary state for UGe2, we now proceed to exam- 
ining electromagnetic properties when the external field 
is applied, regarding UGe2 a conventional type II super- 
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conductor. We note, however, that without H vortices 
with each having unit quantum flux are spontaneously 
created because the internal field ~0.2T far exceeds the 
usual Hci- In this sense the system lacks the complete 
Meissner phase in 77 vs T plane. Let us introduce the 
free energy density fgrad related to the magnetic field H 

fgrad = KiT,a=±^j=y,z\Djri<yf + 7^2 S^^i | £)a;?7o- P (2) 

where Dj = dj — ij^Aj and the unit flux (f)o — We 
have assumed a simple orbital function 0(k) which does 
not break tetragonal symmetry in the basal plane (a, c) 
or (y, z) for simplicity. We reserve its extension to a 
future work when the orbital function turns out to be 
more complex. The essence of the following arguments is 
not altered. Near Tc when H is applied parallel to the z 
axis (the a axis) the above is written as 

hrad = K,{^f + {^fK2{M + ti,Hfx^v+ (3) 
ax 00 

because the magnetic induction B — (O,0,/j,o-ff + M). 
The minimization of the free energy 



00 

readily yields the upper critical field i7°2 as 
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where the transition temperature is redefined as Tj 



(5) 



1ix \/KxK2 



M 



Tc and 



ckQ- When H II y(|| c), by a 



similar way we obtain 



fgrad 



K2{^f + C-^fK,{M^ + (Moff)n^''7+ (6) 
ax 00 



which yields 
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00 



'M2 + (^oi?c'2)' =«o(Tc-T) 
or, approximately near Tc 



(7) 



(8) 



This means that near Tc, 77^2 exhibits a root singularity 
as a function of T with the slope being infinite. This is 
also true for H/ /x{/ /b), namely. 



(9) 



where it shows not only a root singular behavior at Tc, 
but also is scaled with 77^2 by a constant factor. The 
large slope in 77^2 is observed in a certain pressure region 
of UGea by Huxley M 



In the basal (a, c) plane the angular dependence of 
Hc2{S) [6 is the angle from the a axis) is calculated as 



U-(||ncos.^ + (^)^ = L 
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by a simple arithmetics. Since the system is polar axis 
symmetric, namely the magnetization is of vectorial na- 
ture, M > differs from M < 0. Thus when external 
field is reversed from M > to A7 < direction, 77^" 
behaves differently. Namely provided that the magne- 
tization stays its original orientation under the reversed 
field which is 77 < 300G, judging from the magnetization 
curve (Fig. 1 in Ref. 0) 77^*^ is given from eq. (4) 
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for T > Tc- In this situation Tc goes up by 30mK~40mK. 
Then on further increasing field as the magnetization is 
quickly reversed from its original direction and begins de- 
creasing, the Hc2°' curve changes into original 77^2 curve. 
Thus it is expected that the SC state is reentrant by in- 
creasing 77 under a fixed T just above Tc- This Tc rise and 
the associated reentrant phenomenon is deeply rooted in 
the fact that the SC survives under the influence of the 
ferromagnetic state. Upon controlling FM by external 
field we can manipulate and raise Tc- Physically the Tc 
rise occurs because the external field can cancel the spon- 
taneously building internal field due to the ferromagnetic 
polarization which is to lower the hypothetical transition 
temperature and the magnetic induction B can become 
smaller than the internal field A7. 

The possible orbital forms in triplet pairing case al- 
lowed in tetragonal symmetry are given by Volovik and 
Gorkov ||l^ for strong spin-orbit coupling case. Since SC 
appears in the presence of FM, the spin part of the or- 
der parameter is limited to the form: h + ic (a, b and 
c are unit vectors along the three crystal axes). In the 
strong spin orbit coupling case all the one-dimensional 
representations (Ai^, A2u, Bi^ and B2m) are excluded, 
only the pairing function ka(h -f ic) remains as a can- 
didate among the two-dimensional representation E„ in 
their classification. Since the spin orbit coupling is not 
the largest energy scale in the present ferromagnetic su- 
perconductor, the weak spin orbit coupling scheme shown 
next is more relevant classification scheme. 

In the weak spin-orbit case ||l^ , which was the case for 
UPta where another exotic superconducting pairing state 
(either non-unitary bipolar state similar to the present 
material or the triplet planar state) is realized ||l^ there 
are several possible orbital functions allowed coupled to 
the spin function b + ic. Namely, kakckbik"^ — k^) (Ai„), 
kb (A2„), kakch (Bi„) and kb{kl - k^) (B2„) as the 
one-dimensional representations with ki{i — a,b,c) be- 
ing the unit vector for reciprocal space. As for the two- 
dimensional representation E„, the orbital functions of 
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the forms are allowed: Xi{k), Ai(fc) + A2(fc) and Ai(fc) + 
i\2{k) are listed where Ai(A:) = ka and A2(fc) = fcc, or 
Ai(fc) = kaki and A2(fc) = kckl- The bipolar type pair- 
ing state Ai (fc)a-t-?A2(fc)c is excluded. If we take literally 
the orthorhombic crystal symmetry, the strong spin-orbit 
classification does not give rise to a suitable spin function 
of b±ic among the four one-dimensional representations 
such fccC (Ai„), and fcaC, kck (Bi„) and its linear 

combinations. On the other hand, in the weak spin-orbit 
coupling case where there is no two-dimensional repre- 
sentation all one-dimensional presentations have suitable 
orbital functions attached to the spin part h-\-ic., namely, 
kakckb (Ai„), kb (Bi„), k^ (B2„) and /ca (Bsu)- It is in- 
teresting to note that (Ai(A:) -I- iX2(k)){h -\- ic) has not 
only the spin angular moment but also the orbital angu- 
lar moment, both pointing to the ferromagnetic moment 
direction. The allowed states have line node(s) in general 
except for (ka -f ifcc)(b -t- ic) in the two-dimensional rep- 
resentation if the up-spin Fermi surface opens along the 
h axis as predicted by band structure calculation |lC|]. 

We can design several interesting experiments associ- 
ated with the coexistence between FM and SC in the 
zero-field cooling where FM domains are formed, consist- 
ing of M and — M. Each domain has the own SC order 
parameters, either b-|-zc, or b — ic in order to fit its spin 
angular moment with the ferromagnetic moment direc- 
tion. When one of the two domains, say, M with b -I- zc 
percolates throughout a whole system, the macroscopic 
superconducting state establishes. In contrast, under the 
field cooling process whose field strength strong enough 
to make the whole system the single domain, say, M, the 
SC takes place immediately right at T^. 

This reentrant phenomenon mentioned before provides 
direct evidence for the coexistence. The reentrance be- 
havior is indeed observed at 1.35GPa at much higher field 
region (~ a few T) may be different origin from our pre- 
diction. However, it should be kept in mind that the 
superconducting properties in UGe2 show strong sensi- 
tivity of the measured current density |^]. It might be 
related to the domain formation associated with FM. To- 
ward the critical pressure P = 1.7GPa the FM is greatly 
suppressed. There the two long-range orders FM and 
SC truly compete each other while around the optimal 
pressure P ~1.2GPa TpM is much higher than Tc- Thus 
the relationship between FM and SC is changing upon 
increasing P. As a consequence it could be possible to 
change the pairing symmetry from the non-unitary to 
unitary state. This is one of possible explanations of the 
observed reentrance phenomenon at P=1.35GPa. 

We also remark that along the FM domain boundaries 
or the Bloch walls where the magnetization direction ro- 
tates so as to bridge the two oppositely polarized states 
±M, the spontaneously induced spin current flows be- 
cause when transforming h-\- ic state into h — it the SC 
phase also continuously changes. 
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